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Abstract. We determine the dynamical dimer correlation functions of quantum 
dimer models at the Rokhsar-Kivelson point on the bipartite square and cubic lattices 
and the non-bipartite triangular lattice. Based on an algorithmic idea by Henley, 
we simulate a stochastic process of classical dimer configurations in continuous time 
and perform a stochastic analytical continuation to obtain the dynamical correlations 
in momentum space and the frequency domain. This approach allows us to observe 
directly the dispersion relations and the evolution of the spectral intensity within 
the Brillouin zone beyond the single-mode approximation. On the square lattice, we 
confirm analytical predictions related to soft modes close to the wavevectors (n, tt) and 
(it, 0) and further reveal the existence of shadow bands close to the wavevector (0, 0). 
On the cubic lattice the spectrum is also gapless but here only a single soft mode at 
(tt,tt,tt) is found, as predicted by the single mode approximation. The soft mode has 
a quadratic dispersion at very long wavelength, but crosses over to a linear behavior 
very rapidly. We believe this to be the remnant of the linearly dispersing "photon" 
of the Coulomb phase. Finally the triangular lattice is in a fully gapped liquid phase 
where the bottom of the dimer spectrum exhibits a rich structure. At the M point the 
gap is minimal and the spectral response is dominated by a sharp quasiparticle peak. 
On the other hand, at the X point the spectral function is much broader. We sketch 
a possible explanation based on the crossing of the coherent dimer excitations into the 
two-vison continuum. 
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1. Introduction 

Quantum Dimer models (QDM) are important minimal models for the understanding of 
topological order, spin liquid phases and deconfinement. The discovery of a topologically 
ordered resonating valence-bond (RVB) phase of the QDM on the triangular lattice [1] 
has led to an intense exploration of the properties of dimer models in two and three 
dimensions on various lattices [2, 3, 4, 5, 6] and recently also the doping process has 
regained a lot of attention [7, 8, 9]. 

While many ground state properties of QDMs are fairly well understood by now, the 
excitation spectrum has not yet been explored to the same depth. Many basic features of 
the excitation spectrum are only known at the level of the Single Mode Approximation 
(SMA), which gives exact results on the presence of gapless modes. However the actual 
functional form of dimer and other spectral functions are not known in great details. 

Recently, first steps towards a quantitative description of the excitation spectrum 
of the triangular lattice QDM have been undertaken, based on numerical measurements 
of excitation gaps of dimer and vison excitations at [4] and away [10, 11] from the 
Rokhsar-Kivelson (RK) point. The aim of the present study is to go one step beyond 
by determining the full energy resolved dynamical dimer spectral function of QDMs at 
the RK point on three different Bravais lattices in two and three dimensions. 

In two dimensions there are two different classes of phase diagrams of QDMs, those 
for bipartite lattices, such as the square and the honeycomb lattice and those for non- 
bipartite lattices such as the triangular lattice. The bipartite lattices generically allow 
a mapping onto a height model [12] and its subsequent analysis [13] turns out to yield 
critical dimer-dimer correlations and gapless modes at the RK point. The non-bipartite 
lattices instead do not allow such a mapping, and the triangular lattice for example 
has been shown to possess an extended dimer-liquid ground state characterized by 
topological order and a finite gap to excitations [1, 2, 14, 15]. 

In order to highlight the excitation spectra in the different cases, we choose to 
analyze the dynamical dimer correlations on the square and the triangular lattice in two 
dimensions, as well as the three-dimensional simple cubic lattice. 

The plan of our paper is the following. We describe the numerical algorithm, the 
relevant observables and the stochastic analytic continuation method used in Section 2. 
Then we present our main numerical results obtained for the square lattice (Section 3.1) 
and the cubic lattice (Section 3.2), followed by the triangular lattice (Section 3.3) and 
we conclude in Section 4. 

2. Hamiltonian, Observables and the Algorithm 

In the following, we study the simplest versions of QDMs on three different Bravais 
lattices (square, cubic, triangular). The Hamiltonian is defined as follows: 




(1) 



□ 
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where the sum runs over all plaquettes of four sites formed by nearest-neighbor bonds 
only, the Dfup operator flips two parallel dimers on a plaquette, while DFiippabie is a 
diagonal term which is equal to one if the plaquette is flippable, i.e. has two parallel 
dimers, and is zero otherwise. In the following, we focus on the RK point (t = 1, V = 
1) on all lattices, where the equal amplitude superposition of all dimer coverings is a 
ground state [7]. 

In Ref. [16], Henley has proposed a continuous-time Monte Carlo algorithm based 
on local plaquette-flip updates for classical dimers, in which the Monte Carlo time maps 
to the imaginary time axis of the corresponding QDM at the RK point. Therefore by 
calculating Monte Carlo time displaced correlation functions, one equivalently measures 
the imaginary time correlation functions of the quantum dimer model at the RK point. 
This algorithm turned out to be a special case of the more general continuous time 
Diffusion Monte Carlo algorithm [17], where the (quantum) Monte Carlo dynamics 
becomes classical if the ground state wave function is known exactly, as it is the case at 
the RK point. 

This algorithm and its extensions have already been used in the past to extract 
excitation gaps based on a simple fit of the asymptotic long-time behaviour of the 
imaginary time correlation functions, both for dimer- and vison-like excitations on the 
triangular lattice [4, 10, 11]. 

In this work, we now address for the first time the full energy and momentum 
resolved spectral response of dimer excitations on the square, the cubic and the 
triangular lattice by performing a stochastic analytical continuation of the imaginary 
time data collected this way. This procedure allows us to reveal many interesting and 
rich features of the spectral functions, which are beyond hydrodynamic treatments. 

2.1. Observables 

In the following, we focus on the momentum resolved dynamical dimer correlation 
functions. These correlation functions are not only fundamental observables of any 
quantum dimer model, but are also potentially relevant for experiments on spin models 
which map onto quantum dimer models [18, 19, 20]. 

We start by defining the momentum dependent dimer density operator: 



where N denotes the number of sites, Q is a given momentum, r« labels the site 
coordinates and n a (i) counts the number of dimers emanating in the a direction from 
site i . We limit ourselves in this paper to a given dimer direction, here the x axis on 
all lattices. The equal-time dimer structure factor then reads: 



|GS) being the ground state of the Hamiltonian. Finally we define the dynamical 




(2) 



D X (Q) := |iX(Q)|GS>| 2 , 



(3) 
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structure factor as 

V x (Q,u):=--Z lim (GS|4(-Q) D X (Q)\GS) , (4) 

inserting a resolution of the identity gives 

V x (Q,cu) = l(s|4(Q)|GS>| 2 x 6(u - E s ), (5) 

s 

where s labels all the eigenstates with energy E s (Eqs is set to zero). Obviously we 
have D X (Q) = J °° du T> X (Q, uS). A prominent role will also be played by the first energy 
moment of the normalized spectral function 

D _ [fa uV x (Q,u) _ Es E s |(s|^(Q)|GS)| 2 

" fm (Q) : " dm " dm (6) 

2.2. Imaginary time correlations and analytical continuation 

The quantity which is actually measured in the classical Monte Carlo simulations is the 
Laplace transform of (4), resp. (5): 

/>oo 

T>x(Q,t) = / dcjexp[-wr] V x (Q,u) (7) 
= |{s|^.(Q)|GS>| 2 x exp[-£?.r]. (8) 

s 

The difficult part is then to recover the spectral weights |(s|D x (Q)|GS)| 2 and the 
excitation energies E s , based on Monte Carlo data which inevitably contain error bars. 
This is a well-known example of an ill-conditioned problem and originates from the fact 
that time-displaced observables are obtained in imaginary time. In the Monte Carlo 
community, this problem has been addressed to some extent by enforcing a Maximum 
Entropy criterion to get dynamical spectra [21]. The presence of an entropy term in 
the minimization procedure allows one to obtain the most probable spectrum given the 
Monte Carlo data and has been shown to be reasonably accurate in many situations. In 
this work, we use a sophisticated formulation called Stochastic Analytic Continuation 
(SAC) initiated by Sandvik [22] and formalized in detail by Beach [23], where instead 
of obtaining a single spectrum, many "acceptable" spectra are averaged. There is 
some freedom left in the averaging procedure, but we emphasize that our findings and 
conclusions are robust with respect to this freedom. For instance, although sum-rules 
can be enforced in the Maximum Entropy procedure, we refrain to do this, and instead 
we systematically check that the first moment of the spectral distribution obtained by 
SAC is identical (within the error bars) to the first moment obtained directly by the 
Monte Carlo simulations by virtue of the relation o;^(Q) = — 4-T> x (Q, r)/D x (Q). Via 
this procedure we obtain an independent assessment of the quality of the stochastic 
analytical continuation. 
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2.3. Single mode approximation (SMA) 

On the analytical side, the single mode approximation (SMA) has proven to be a useful 
tool for the analysis of some aspects of the excitation spectrum of QDMs. Based on its 
definition 

d x (n , l (GS\[D x (-Q),[H,D x (Q)]]\GS) f . 

-sma(Q) •= 2 dM > (9) 

and the fact that this quantity constitutes an upper bound for the dimer gap at 
wavevector Q, the absence of a gap can simply be shown by demonstrating that either 
the numerator vanishes or the denominator diverg function of system size. In 

the seminal paper by Rokhsar and Kivelson [7] they indeed used the fact that the 
numerator vanishes at (tt,tt) (D X (Q) with Q = (tt, tt) is a conserved quantity) to prove 
the gapless nature of the square lattice dimer spectrum at (it, it). Similarly, one can 
prove the absence of a dimer gap at (tt, tt, tt) for the simple cubic lattice at the RK 
point [5, 6]. Furthermore, for square lattice again, Moessner and Sondhi [6] showed that 
the logarithmically divergent structure factor D X (Q) with Q = (tt, 0) leads via the SMA 
to another soft mode, which they called piOn. 

Another important property of the SMA is that it is equivalent to the exact first 
frequency moment of the dynamical correlation function which is used to construct the 
trial excited state. In our case this is precisely the dimer density operator at wavevector 
Q and a given orientation. Then 

1 f°° 

^sma(Q) = D -^ ) J Q duj oo V x (Q,u) = w£(Q). (10) 

Historically the SMA bears its name from the approximation 

V S X MA (Q, U ) » D X (Q) x 5(uj - wg^CQ)), (11) 

which is very sucessful for a large number of systems, most prominently for 4 He as shown 
by Feynman [24], but also for quantum Hall systems [25] and Haldane spin chains [26]. 
As an important result, we will show below that, for the QDMs, the SMA gives in general 
a rather poor account of the features contained in the dynamical dimer spectrum, with 
the exception of the hydrodynamic regime in the close neighborhood of the (tt, tt) point 
on the square lattice and the (tt, tt, tt) point for the simple cubic lattice. 



3. Numerical results for the dynamical dimer correlation functions 

In this section we present our numerical spectra for the various lattices. The Monte 
Carlo data has been obtained on systems of linear size L = 64 for the square lattice, 
L = 32 for the cubic lattice and L = 36 for the triangular lattice. For the spectral 
function measurements we used a uniform grid of At = 0.1 up to r max = 80 for the 
square and the cubic lattices. For the triangular lattice the required r max value was 
smaller due to the finite gap. For precise measurements of the r derivative at the origin 
finer grids in At have been used. 
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Figure 1. Square Lattice. Upper panel: dynamical dimer spectral functions 
T> x (Q,u) plotted along the Brillouin zone path (0, 0)-(n, 0)-(tt, 7t)-(0, 0) on a system 
of linear extension L — 64. The first moment cj fn ^(Q) of the distribution - equivalent 
to the SMA prediction - is shown by the white filled circles. Lower panel: equal-time 
dimer structure factor D X (Q) along the same path. The peak at (it, 0) is logarithmically 
divergent with system size. 



Our single plaquette flip dynamics conserves the winding number on the square 
lattice, the fluxes on the cubic lattice and the topological sectors on the triangular 
lattice. We chose then to work in the zero winding number / zero flux sector for the 
square and cubic lattices respectively, while we worked in the trivial topological sector 
for the triangular lattice. For the square and the cubic lattice this choice might induce 
a slight bias towards the physics at v/t < 1, since there the zero winding number / flux 
sector is preferred. On the triangular lattice however the topological degeneracy insures 
that the results in the different topological sectors are identical. 

3.1. Square Lattice 

The RK point of the square lattice is a special point, since it is a critical point separating 
a staggered dimer phase at larger values of v/t > 1 from a symmetry broken phase at 
v/t < 1, whose detailed structure is still a topic of ongoing research [27, 28, 29, 30]. 
Since the seminal work of Rokhsar and Kivelson [7], it is known that the excitation 
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Figure 2. Square Lattice. Left panel: low-energy zoom ofD x (Q,ui) into the (n,0) 
region. Right panel: similar zoom into the (tt, tt) region. The dashed white lines with 
the filled circles denotes the measured u)^(Q) for both panels. The dash-dotted line 
in the left panel is the analytical prediction (15) for the maximum intensity at low 
energy, while in the right panel the dash-dotted line shows the analytical prediction 
(12) for the resonon dispersion. 



spectrum is gapless and has a quadratically dispersing dimer-wave mode becoming soft 
at (7r,7r). This SMA dispersion is given by [7, 31]: 

^resonon(k) = -k 2 + higher order terms, (12) 

o 

where Q = k + (tt, tt), and k is the modulus of k. It was realized that (tt, tt) is not the 
only gapless point in the Brillouin zone, when Ref. [6] characterized a novel excitation 
termed piOn which is gapless at (tt, 0) and (0,tt), and which has its root in the height 
model description [13]. However despite this knowledge of the " long- wavelength" , low- 
energy behaviour, a full characterization of the dimer spectral response is lacking. 

We present in the upper panel of Fig. 1 our numerical results for the dynamical 
dimer structure factor V x (Q,u) at the RK point on the square lattice obtained by the 
numerical simulations outlined in section 2. The color highlights the intensity at a given 
Q and on, while the line with the white circles denotes the first moment oo^(Q). In the 
lower panel the equal-time structure factor L\.(Q) is shown, which agrees nicely with 
the analytical results of [6, 31] (not shown). Note that since the parity of the number 
of dimers crossing any vertical line is conserved by the dynamics, [D x (q, tt), Hqdm] = 0, 
which implies that the spectral weight vanishes along the Q = (q, tt) line, including the 
(tt, tt) point. The non-continuity of the equal-time structure factor upon approaching 
(tt, tt) from different directions reflects the pinch (bow-tie) singularity in momentum 
space [32], and is a generic feature of two-dimensional bipartite Rokhsar-Kivelson 
points, giving rise to the critical real space dimer correlations. On the other hand 
the logarithmic divergence of D x (tt, 0) is specific to some models, since a lower stiffness 
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K in the low-energy action can lead to a non-divergent equal-time response, as is the 
case for example in a generalized classical model [32]. As expected the SMA prediction 
matches the full spectral function very well in the close vicinity of the (jr, tt) point. 
However in the rest of the Brillouin zone the intensity is by no means concentrated near 
the SMA result. We can however identify traces of intensity above the SMA curve, 
loosely following it, for example in the neighborhood of T — (0, 0). On the other hand 
there is also some intensity at low frequencies and approaching zero energy as Q —>■ T. 
Further investigation is needed to clarify whether this low energy signal is a "shadow 
band" feature of the (tt, 0) or the (tt, tt) structure. 

We now focus on the low-energy features around the (vr,0) and the (tt, tt) region, 
which are presented in a zoomed view in the two panels of Fig. 2. First we briefly 
discuss the behavior in the (tt, tt) region shown in the right panel. As stated before the 
SMA predicts a quadratically dispersing mode in the vicinity of this point. Indeed our 
numerical results follow the refined analytical prediction (12) presented by the dashed 
line very closely and the whole spectral function is actually dominated by the weight 
at this energy. It is difficult based on our numerical data for the spectral functions to 
address the linewidth of the mode, but based on the following argument we expect a 
finite lifetime for this quadratic mode. As a general remark for a quadratic dispersion, 
it can easily be understood that for any Q, the gap should vanish in the thermodynamic 
limit. Indeed, by combining many low- lying excitations, such that the total momentum 
is Q, an arbitrary small energy state can be obtained. An equivalent statement is that 
a given Q excitation could decay into many lower-energy excitations (when there is no 
symmetry preventing this process). Therefore, if we observe a quadratic behavior over 
a large region, this is because the spectral weight on the lower energies states is tiny, 
due to density of states and matrix element effects. 

The long-wavelength, long-time behavior of the (tt, 0) feature has been discussed in 
Ref. [6]. Based on the analysis of the height model representation the asymptotic form 
of the imaginary time correlation is known in the vicinity of (tt, 0) 

/(k,r)~exp(-AVF), (13) 

where now Q = k+(7r, 0), and k is the modulus of k. This expression can be transformed 
by an inverse Laplace transform giving: 

n^) = k -^M- d4) 



2yW/ 2 

The expression (13) is supposed to be valid only at long times, as can be seen from the 
fact that the first moment of -F(k, uj) does not exist due to the infinite r derivative of 
/(k, r) at r = 0, while the true dynamical dimer structure factor has a bounded first 
moment for all k. Even though the first moment does not exist, it is possible to track 
the maximum of F(k, u) as a function of k, giving the simple expression: 

uj pi0n (k) = h 2 . (15) 
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Figure 3. Cubic Lattice. Upper panel: dynamical dimer spectral functions 
T> X (Q, lj) plotted along the Brillouin zone path T — (0, 0, 0) — > (tt, 0, 0) — > (tt, n, tt) — > T 
on a system of linear extension L = 32. The first moment ui^* (Q) of the distribution 
- equivalent to the SMA prediction - is shown by the white filled circles. Lower panel: 
equal-time dimer structure factor D X (Q) along the same path. Note that D X (Q) has 
no divergent feature on the cubic lattice. 

Note that the location of the maximum is also close to the onset of intensity in the 
formula (14). We plot this simple result by a dash-dotted line in the left panel of Fig. 2, 
and find very a nice agreement between this parameter-free prediction and the actual 
numerical results. The dimer excitation spectrum has definitely no quasiparticle peak 
in the vicinity of the (tt, 0) point, but instead a characteristic broad continuum. 

3.2. Cubic Lattice 

We now turn to the three-dimensional (3d) cubic lattice. Being also bipartite, a well 
defined mapping to a long wavelength action exists [5], and this model is expected 
to be somewhat similar to the two-dimensional square lattice case. In particular the 
SMA argument predicts a quadratic dispersion close to the (n, tt, it) point [6]. However, 
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Q Q 

Figure 4. Cubic Lattice. Left panel: low-energy zoom u> S [0,0.3] of T> x (Q,lu) 
along the same path in the Brillouin zone as in Figure 3. Right panel: Supplementary 
low energy data for ^(Q,^) along the path (ir,ir,0) — > (tt,it, tt). The dashed white 
lines with the filled circles denotes the measured iv®* (Q) for both panels, and is in good 
agreement with the expected k 2 dispersion close to (7r,7r, tt). Note that the maximum 
intensity of 2? x (Q,oj) rapidly deviates from the expected fc 2 -law and becomes linear, 
revealing the remnant of the linearly dispersing "photon" mode of the J7(l)-liquid 
Coulomb phase. 

contrary to the square lattice case, there is no other peculiar point in the Brillouin zone 
(such as the 2d "piOn") so that the (7r,7r,7r) point is supposed to be the only gapless 
excitation. Also, in the 3d case, there is a so-called Coulomb phase with dipolar dimer 
correlations [5, 6, 27]. Moreover, based on an effective low energy field-theory [5, 6], 
the cubic dimer model can be described: (a) for v c /t < v/t < 1, by a standard U(l) 
Maxwell electromagnetic theory with a linearly dispersing transverse photon; (b) at the 
RK point, v = t, the photon dispersion becomes quadratic. 

Fig. 3 shows our results both for the dimer spectral functions and equal-time 
structure factors along a specific path through the Brillouin zone. Except at (it, tt, tt) 
where it vanishes by symmetry, the equal-time structure factor D X (Q) is rather smooth 
and has no divergence. The non-continuity of D x (Cl) upon approaching (tt, tt, tt) from 
different directions is the momentum space signature of the dipolar real-space dimer 
correlation functions. 

The first moment of the spectral functions, which corresponds to the SMA result, 
exhibits a wide dispersion reaching its maximum of ~ 2.45 at the T point, with a 
quadratic behavior aproaching zero energy close to (tt, tt, tt), similar to the square lattice. 
The visible intensity is however concentrated mostly at very low energy, in contrast to 
the square lattice. By looking more closely at the low-energy features shown in Fig. 4, 
it can be seen that the spectral function is not simply following the SMA energy, even 
in the neighorhood of (tt, tt, tt), clearly at variance with the results for the square lattice, 
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where the intensity was well described by the SMA prediction near (vr,7r), see the right 
panel of Fig. 2. Instead we see a linearly dispersing feature with an approximate velocity 
c ~ 0.04(1) in all directions explored around (jr, tt, tt). Of course, such a linear behavior 
cannot exist arbitrary close to (tt, it, tt), because the SMA correctly predicts a quadratic 
dispersion, but nothing prevents such a linear mode with a tiny velocity further away 
from this point, i.e. for energies above a certain crossover scale. This unexpected feature 
probably has its origin in the truly linearly dispersing photon mode expected at smaller 
v/t adjacent to the RK point, and whose velocity is supposed to grow as a function of 
(1 — v/t). Our observation calls for an investigation of the cubic QDM using Green's 
function Monte Carlo techniques in order to explore the dispersion of the photon deep 
inside the Coulomb phase and when approaching the RK point from smaller v/t < 1 
values. 

3. 3. Triangular Lattice 

In contrast to the previous two bipartite lattices, the triangular lattice QDM at the RK 
point is known to be in a gapped Z 2 -liquid phase sustaining topological order [1], which 
is is stable in a finite region of parameter space for 0.8 <v/t< 1. The dimer excitations 
possess a finite gap in the whole Brillouin zone [4]. 

Our numerical data shown in Fig. 5 clearly reveal the presence of a finite gap. 
Even more striking is the fact that the first spectral moment lu^(Q) is nowhere in the 
Brillouin zone close to the the intensity carrying low energy excitations, and even the 
location of the minimum of the SMA predication (K point) and the actual lowest dimer 
excitation (M point) do not coincide. This is in contrast to the square and cubic lattice, 
where the SMA gave a faithful description of the dimer exciation energies at the (it, tt) 
and the (tt, tt, tt) point respectively. The equal-time structure factor is smooth on the 
triangular lattice, due to the short-ranged nature of the dimer correlations in real space. 

Looking more closely at the low energy part of the spectral functions in Fig. 6 one 
recognizes that the intensity carrying modes at the onset of the dimer spectrum are much 
less dispersive than the SMA prediction suggests. Furthermore the broadening can vary 
through the Brillouin zone. For instance, the spectrum is very broad at the X point, 
while it is much sharper close to the M point, see Fig. 7. These observations strongly 
suggest that the dimer waves excited by the operator (2) are not the true "elementary" 
excitations of this dimer model. 

Indeed the elementary excitations are known to consist of so-called visons which 
are non-local in terms of dimers and for which the dispersion has been computed [4]. 
Real dimer excitations can be understood as 2- or more vison excitations. Here, the 
lowest-energy dimer excitation (at the M point) is formed by a bound-state of two 
visons [4]. Based on the single vison dispersion relation [4], we have plotted in Fig. 6 
the bottom of the 2- vison continuum. By comparing to the spectral weight of the dimer 
correlations, we observe that: (a) close to K and T points, the spectral weight nicely 
follows the 2-vison continuum; (b) the lowest energy state at M is slightly below the 
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Figure 5. Triangular Lattice. Upper panel: dynamical dimer spectral functions 
2? X (Q, ui) plotted along the Brillouin zone path: T^X^K^M— > T for a system 
of linear extension L = 36. The points in the Brillouin zone are labeled in the inset. 
The first moment u>F* (Q) of the distribution - equivalent to the SMA prediction - 
is shown by the white filled circles. Lower panel: equal-time dimer structure factor 
D X (Q) along the same path. 

two-vison gap, which is consistent with the interpretation as a two-vison bound-state; 
(c) at the X point, the spectral weight is very broad, we cannot reveal any well-formed 
excitation so that a dimer excitation probably disintegrates in the 2-vison continuum, 
a point of view which has also been raised in Refs. [4, 11]. These observations then 
suggest a scenario where the coherent (quasiparticle-like) dimer excitations close to the 
M point cross into the continuum going either to the K or the V point, giving rise to 
damped, incoherent excitations for example at the X point. It is also possible that there 
are termination points similar to magnetic systems or liquid Helium 4 [33]. 

A natural extension of this work would be to investigate the spectral functions of 
the visons directly, first at the RK point [4], and then also towards the phase transition 
from the topologically ordered phase into the symmetry broken a/12 x V4~2 phase [11]. 
It will be interesting to confirm that the visons have a sizable quasiparticle weight in the 
topologically ordered phase, and therefore can really be considered as the elementary 
excitations inside this phase. Also the currently unexplored nature of the excitations 
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Figure 6. Triangular Lattice. Low-energy part of the dynamical dimer spectral 
function along the same path in the Brillouin zone as in Figure 5. The onset of the 
two-vison continuum is indicated by the white stars. The single vison energies used 
here have been obtained in Ref. [4]. 




Figure 7. Triangular Lattice. Comparison of Classical Maximum-Entropy and 
Stochastic analytical continuation spectra for the dimer spectral functions at the X and 
the M points. Both methods for analytic continuation yield very similar spectra. The 
vertical green dash-dotted lines indicate the dimer gap energies extracted in Ref. [10]. 
Insets: Monte Carlo imaginary time data used for the analytical continuation. The 
dash-dotted green line indicates the fit with a single exponential decay obtained in [10]. 
It is obvious that the spectral function at the X-point has no sizable quasiparticle 
peak, while the spectral function at the M-point seems to show a resolution limited 
peak at the lowest energy, providing evidence for a coherent dimer excitation at this 
momentum. 
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above the ^/l2 x y/Y2 ordered groundstate (dimers vs visons) could be clarified this way. 

4. Conclusions and perspectives 

Quantum dimer models at the Rokhsar-Kivelson point can be efficiently simulated with 
a classical Monte-Carlo algorithm and in particular, imaginary time-displaced dimer 
correlations can be obtained. By performing a stochastic analytic continuation to real 
frequencies, we have been able to compute quantum dimer spectral functions for a 
variety of lattices. 

For bipartite lattices (2d square and 3d cubic), we have confirmed that these 
models are in a gapless phase, with a rather well-defined resonon excitation dispersing 
quadratically. The 2d case also exhibits other gapless excitations close to (tt, 0) (piOns), 
that are not described by a simple dispersion but rather a broad spectrum in full 
agreement with analytic calculations. In the 3d case, on top of the quadratic dispersing 
mode, we have shown evidence that a linear-dispersing excitation might play a role at 
shorter wavelengths and might be related to the well-known photon-like excitation that 
appears away from the RK point. 

For non-bipartite lattices (2d triangular), these models are in a Z 2 -liquid phase 
with a finite-gap for dimer excitations. By analyzing the low-energy structures, we have 
provided arguments that, for different points in the Brillouin zone, the lowest energy 
dimer excitations can be viewed either as a coherent two-vison bound-state, or on the 
contrary can decay into the two-vison continuum. 

Similar future studies could investigate further open questions, just to cite a few : 

(i) analytic continuation of Green's function Monte Carlo data away from the RK point 
(the possibility to approximately compute time-displaced correlations has already 
been used in Ref. [10]), observation of the linearly dispersing "photon" of U(l) 
gauge theories on 3d bipartite lattices [5, 6, 27]. 

(ii) spectral functions and coherence of the visons on the triangular lattice. 

(iii) The stiffness appearing in the low-energy approach can be varied by changing the 
model. In Ref. [34], it has been shown that, by adding an aligning term, the stiffness 
can be increased continuously from 1/2 to 4, before a phase transition occurs. In 
Ref. [32] , by including longer-range bipartite dimers, the stiffness changes smoothly 
from 1/2 to 1/9, while allowing for non-bipartite dimers immediately opens a gap. 
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